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Abstract 

With the aim of resolving theoretical issues associated with the fourth root 
prescription for dynamical staggered fermions in Lattice QCD simulations, we 
consider the problem of finding a viable lattice Dirac operator D such that 
(detD staggered) 1 ^ = detD. Working in the flavour field representation we show 
that in the free field case there is a simple and natural candidate D satisfying 
this relation, and we show that it has acceptable locality behaviour: exponen- 
tially local with localisation range vanishing ~ yja/m for lattice spacing a — > 0. 
Prospects for the interacting case are also discussed, although we do not solve 
this case here. 
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1 Introduction 



The development in recent years of an improved staggered fermion formulation [1] has 
made unquenched numerical Lattice QCD simulations possible at realistically small 
quark masses. The resulting impressive agreement between the calculated parame- 
ters of QCD phenomenology and their experimental values [2, 3] (along with predic- 
tions for quantities not yet measured experimentally [4]) indicates that the long-time 
dream of being able to do high-precision Lattice QCD calculations is now becoming 
reality. However, the advantageous properties of staggered fermions for numerical 
implementation are currently offset by unresolved problematic issues at the concep- 
tual/theoretical level. In particular, there is concern [5, 6, 7] about the use of the 
fourth root of the staggered fermion determinant to represent the fermion determinant 
of a single dynamical (sea) quark. A number of works have appeared recently address- 
ing this concern via theoretical considerations [8, 9, 10, 11, 12, 13, 14, 15], derivations 
of predictions that can be used to test its viability [16], and various numerical investi- 
gations [17, 18, 19]. The present paper is intended as another theoretical contribution 
in this direction. 

A staggered fermion is a lattice formulation of four continuum fermion flavours, 
nowadays called "tastes" (to distinguish them from the actual quark flavours). The 
fermion determinant for a single quark flavour in this framework is represented by 
a rooted determinant det(D 'staggered) 1 ^ 4 ■ While this formally goes over to the deter- 
minant for a single quark flavour in the continuum limit, the concern regarding this 
prescription is that it does not fit in an obvious way into the framework of local lattice 
field theory at non-zero lattice spacing. The lattice model might therefore not be in 
the right universality class to reproduce QCD. This raises the question of whether 
the dynamical staggered fermion formulation is a first principles approach to QCD or 
simply a phenomenological model which describes QCD very well in a certain regime. 

One way to establish that the universality class is the right one would be to show 
that there is a viable (and in particular, local) single-flavour lattice Dirac operator D 
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such that [5] 



(detD staggered ) 1/4 = detD 



(1.1) 



This would imply equivalence between the dynamical staggered fermion formulation 
and the manifestly local formulation with sea quarks described by D. (To avoid 
unitarity issues, D, with suitably adjusted bare mass, should then also be used as 
the Dirac operator for the valence quarks.) We will refer to (1.1) as the Staggered 
Determinant Relation (SDR) in the following. 

The most direct attempt at a solution to the SDR is simply to take D = (D 'staggered) 1 ^ ■ 
This is essentially the approach taken by Jansen and collaborators in Ref.[9]. More 
precisely, they considered the operator 



for which detM = det(D staggered ) l/2 since (D staggered y D staggered couples lattice sites 
by even-to-even and odd-to-odd. Thus M is a candidate operator for D in the case 
where 1/4 — > 1/2 in (1.1), i.e. the case of two degenerate quark flavours. However, 
this operator was found to have unacceptable locality behaviour: it is exponentially 
local (for bare mass m > 0), but the localisation range is ~ mT 1 and thus fails to 
vanish in the limit of vanishing lattice spacing, a — > [9]. 

This negative result of Ref. [9] is unsurprising, since the operator (1.2) does not 
take account of the staggered fermion taste structure. The staggered fermion action 
can be viewed as consisting of naive fermion actions for four fermion species (the 
tastes), together with terms that couple these, with the latter formally vanishing for 
a — > 0. This suggests that, in attempting to find a local, single-flavour lattice Dirac 
operator satisfying the SDR, one should consider operators of the form 



where yl is the (massless) naive lattice Dirac operator, and W is a term which formally 
vanishes for a — > and whose role is to take account of the taste-mixing in D staggered . 
In this paper, working in the flavour field representation of staggered fermions [20], we 




(1.2) 



D =yf +W + m 



(1.3) 
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show that, in the free field case, there is a natural candidate D of the form (1.3) which 
satisfies the SDR. The W operator, although formally vanishing for a — > 0, turns out 
to involve a square root, so its locality status is not immediately clear. We show, 
however, that the operator does have acceptable locality behaviour: exponentially 
local, with localisation range vanishing ~ y 'a/m for a — > 0. 

Our operator can be gauged (i.e. coupled to the link variables) in a variety of 
ways. However, for reasons which we will discuss later, it is most unlikely that a 
gauging of this operator exists such that the SDR continues to hold in the interacting 
case. Our operator should therefore be regarded as a prototype, i.e. a first step on a 
path to constructing more sophisticated operators which have a chance of satisfying 
the SDR in the full interacting case. 

Of course, there is no a priori guarantee that a viable lattice Dirac operator 
satisfying the SDR actually exists in the interacting case, so other approaches should 
also be considered. One possibility is the following: If there is a single-flavour D 
such that the effect of including the determinant ratio det(D sta ggered) 1 ^ / detD in the 
Lattice QCD functional integral is simply to renormalise the bare coupling constant 
(just as dynamical heavy quarks do [21]), then representing the sea quark determinant 
by det{D staggered) 1 ^ is equivalent to representing it by detD together with a shift in 
the bare coupling. Since the latter description is manifestly local, this is another 
way in which the locality issue could be positively resolved. The prospects for this, 
and the properties that such a D would be expected to have, are also discussed in 
some detail in this paper. The problem of finding such a D is seen to be essentially 
equivalent to the problem of finding a solution to a generalised version of the SDR. 

The paper is organised as follows. After a general discussion of the problem of 
finding viable solutions to the SDR, we arrive at our free field candidate D in §2. 
In §3 we prove that this operator has the good locality behaviour mentioned earlier. 
Our argument is entirely analytic and the techniques are of a generally applicable 
nature; we also apply them to give a new derivation of the negative locality result for 
the operator considered previously by Jansen and collaborators [9] . (Their argument 
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in the free field case had a numerical as well as analytic component.) We conclude in 
§4 with a discussion of the issues and prospects for the interacting case. 

2 Considerations for finding a candidate D 

For concreteness we specialise to 4 spacetime dimensions in this section (everything 
generalises straightforwardly to arbitrary even dimensions). The usual staggered 
fermion action, obtained via spin-diagonalisation of the the naive action [22], is 
Sstaggered = a 4 J2x "4> (x) D st ip (x) where the staggered Dirac operator is given by 

D st = ^-V^ + m (2.1) 

(X 

with rj^(x) = (_ i)( :E i+---+ ;c M-i)/ a a nd = |(V+ + V~) the symmetrised gauge- 
covariant difference operator. The usual flavour (i.e. taste) identification comes about 
by considering the free field propagator: it has 4 poles, and the momentum space 
Brillouin zone is divided into 4 subregions, each containing a single pole, with the 
momenta in each of these subregions being interpreted as the momenta of different 
fermion tastes. An alternative, and conceptually more appealing way to identify 
the tastes is provided by the flavour (taste) field representation of the staggered 
fermion action derived in [20]. In this representation the taste fields are manifest 
from the beginning in the fermion action. The taste fields live on the blocked lattice 
(spacing=2a), whereas the lattice paths and link variables which specify the gauging 
of the action are those of the original lattice. The action in general gauge background 
does not have a simple expression in this setting though, making it more difficult to 
work with in practise. However, in the free field case the action does have a simple 
expression. Denoting the staggered Dirac operator in the taste field representation 
by D stt , it can be written in the free field case as [20] 2 

Dlu e = ® l)^V p + *( 75 ® A, + m (2.2) 

2 The free field version of the staggered Dirac operator in the taste representation also arises from 
a first principles approach to constructing the Dirac operator on the lattice [23] . 
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where now the (free field) difference operators are on the blocked lattice; A u = V~ — 
(so that A = J2 U is (2a) 2 times the blocked lattice Laplace operator), and {r^} 
is a hermitian representation of the Dirac 7-algebra on taste space C 4 . 3 

The importance of taking account of taste structure when attempting to find a 
solution D to the SDR can now be seen in the free field case as follows. The 7-matrix 
representation {T v } in (2.2) can be chosen such that the diagonal matrix elements in 
each of the P^s all vanish. Then the taste-mixing terms in the free field lagrangian 
correspond to the terms with P^s in the free field Dirac operator (2.2). Therefore, 
if the taste-mixing terms are "turned off" the free field Dirac operator reduces to 
y <S> 1 + m, where y = 7 M ^V M is the (massless) free field naive Dirac operator and 
1 is the identity matrix on taste space. Consequently, 

(detD{n 1/2 - det{^ ^y^J , (detD^) 1 ^ - dettf + m) (2.3) 
or, alternatively, 



(detDir^ - det N& + rnW + m) 
V st ! \ v /(y + m)t(y + mj 

(detD f s r) 1/4 -> det(y/Qt + myQ# + m) ) . (2.4) 



In the former case the fractional powers of detDg[ ee become determinants of ultra- 
local lattice Dirac operators, while in the latter case they become determinants of 
operators which cannot be expected to have good locality properties. If we now 
imagine turning back on the taste-mixing terms, there is reason to hope that there 
will be corresponding deformations of ^V+ m y/* m ^ or y + m into some two-taste 
lattice Dirac operator D or single-taste D, respectively, which continues to have good 
locality behaviour, such that (detD s s r t ee ) 1 ' 2 = detD and (detD s J t ee ) l / A = detD. On 
the other hand, if a solution D to the SDR, or a solution D to the version of the 
SDR with with fractional power 1/2 of the staggered fermion determinant, has been 

constructed "blindly" without taking account of the taste structure of the staggered 
3 The part 175 (giT" in (2.2) is usually written as 75(g)T 5 r ly where {t v } is a hermitian representation 
of the Dirac 7-algebra on taste space. Note that T u = —it^t" defines another (equivalent) hermitian 
representation of this algebra. 



fermion formulation, it can happen that when taste-mixing terms are turned off in 
the free field case the scenario (2.4) arises; then it is to be expected that the D or D 
have bad locality behaviour. In fact this is essentially the situation for the solution 
D = M considered in [9], and the negative locality result found there is therefore 
unsurprising. However, the possibility (2.3) gives hope of doing better than this, at 
least in the free field case. 4 

In light of (2.3), when attempting to find a viable D in the free field case it is 
natural to consider Dirac operators of the form 

c ^4 v " + i w/+m (2 - 5) 

on the blocked lattice, where the purpose of -^W is to take account of the taste- mixing 
terms in the staggered Dirac operator. In particular, W should formally vanish ~ a 2 
for a — > 0, and should lift the species doubling of the naive Dirac operator. In other 
words, -^W is to be a Wilson-type term. 

A feature of the free field staggered Dirac operator (2.2) in the taste field rep- 
resentation, which is very useful in this context, is that (-D^ ee )^ D{[ t ee is trivial in 
spinor®flavour space: 

(D5n*D£° = ( -V 2 + E(!^) 2 + (2am) 2 ) (1 ® 1) (2.6) 

On the other hand, for a free field operator of the form (2.5) we have 

D^D = -^(-V 2 + (W + 2am) 2 )l (2.7) 

trivial in spinor space. Comparing (2.6) and (2.7), and noting that detD = det(D^ ' D) 1 ^ 2 
(assuming + m > 0) and det(D st ) = det(D\ t D st ) 1 / 2 (assuming m > 0), we im- 
mediately see that a sufficient criteria for the desired determinant relation detD = 
{detD f J ee y/ 4 to be satisfied is 

(W + 2am) 2 = ^(|A^) 2 + (2am) 2 . (2.8) 

V 

4 The interacting case is more difficult, since the taste field representation of the staggered Dirac 
operator is not given simply by some gauging of the V^'s and A^'s in (2.2) but has a more compli- 
cated structure [20] . 
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This has the solution 



W = /(2am)2 + ^(|A^ 



— 2am , 



(2.9) 



which clearly has the required properties for -^W to be Wilson-type term (i.e. W 
lifts species doubling and formally vanishes ~ a 2 for a — > 0). Thus we have arrived 
at a free field solution D to the SDR (1.1). Substituting (2.9) into (2.5) we get the 
expression 



Note that turning off the taste-mixing terms in the free field staggered fermion action, 
which, as pointed out previously, corresponds to putting Y v — > in (2.2), has the same 



a realisation of the scenario (2.3). However, because of the square root in (2.10), it 
is not immediately clear that good locality behaviour of D, anticipated in our earlier 
discussion, is realised. In fact this square root operator has some similarity with the 
free field square root operator considered by Jansen and collaborators in [9], which 
turned out to have unacceptable locality behaviour. Nevertheless, we show in the 
next section that our operator does have good locality behaviour. The reason why 
its behaviour is different from the operator in [9] is a bit subtle, and to elucidate this 
we also provide in the next section a new derivation of the negative locality result of 
[9] which reveals the origin of the different behaviours. 

We remark that other, ultra-local solutions to the SDR exist in the free field case. 



Using -V 2 + (±A,) 2 = A u (2.6) reduces to [D{T)^{T = (t^A + m 2 )(l ® 1) 



and it follows that the free field SDR is satisfied, e.g., by D = + m 2 ) 2 acting 

on scalar Grassmann fields on the lattice. Other examples of ultra-local solutions are 
easily constructed. However, these are unattractive options since they do not have 
the form of a lattice Dirac operator. 




(2.10) 



effect as putting A„ — > 0. By (2.10) this gives D -±yl + m (for m > 0); thus we have 



8 



3 Free field locality result 



In this section we work in arbitrary spacetime dimension d and show that the free 
field operator 



^(amy + ^Al (3.1) 
on lattice with spacing a is exponentially local with localisation range vanishing ~ 



^Ja/m for a — > 0; then D in (2.10) obviously has this same locality behaviour on the 
blocked lattice. The argument proceeds in several steps. First, we specialise to d— 1 
dimension and write 

(am) 2 + A 2 (x,y) = -T^ J (am) 2 + A(p) 2 e ip{x - y)/a 

= _L / ^ h am )2 + {2 -(z + Z-^)) 2 Z \ X - y \l a 

2maJ\ z \=i z vv ; v v 
dz 



i T v z ~ 2 (( am ) 2 "* 2 + (2 ~ 1)4 ) - 2 '^ l/a (3 - 2) 



27ria /|z|= 

(the integral is counter-clockwise around the unit circle in the complex plane and we 
have set z = e ±tp with "+" if x — y > and "— " if x — y < ; for x = y either choice 
can be used.) The square root z \— > \fz is holomorphic after making a cut in C ; we 
choose the cut to be the half-line of negative real numbers R_ . Then, by the residue 
theorem, the circle around which the integral in (3.2) is performed can be shrunk to 
a closed loop around the region containing the z's for which f(z) G R_ , where 

f(z) = z- 2 ((am) 2 z 2 + (z-l) 4 ) (3.3) 

is the function inside the square root in (3.2), since outside this region (and away 



from z — 0) \J f{z) is holomorphic. The excluded z's are are found as follows: 

f(z) = -A A G R + 
& (z- l) 4 + ((am) 2 + X)z 2 = 

^ [(z - l) 2 + i^J(am) 2 + A z) {{z - I) 2 - (am) 2 + A z) = (3.4) 

For given A G R + there are 4 solutions; we are only interested in the ones with \z\ < 1 
and these are z — z± (s = ^f(am) 2 + \) where 

ts i 

z± = 1 ± — - y/syJM - s/4 s G] 0, oo[ (3.5) 
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Figure 1: The "exclusion curves" C±. The locations of the endpoints near 1 depend 
on am and converge to 1 for a — > 0. 

Thus the z's for which /(z) G R and |z| < 1 form curves inside the unit circle in C, 
parameterised by (3.5) with s G [am, oo[ . It is useful to re-parameterise these curves 
as follows. We introduce 



t = 1 



s/2 



(3.6) 



note that this is a strictly decreasing function of s with t — 1 for s = and t — > for 
s — > oo . After a little calculation (3.5) can be re-expressed in terms of t as 



z± (t)=tTi(l-t)^^- i te]0,l[ (3.7) 

The z's for which /(z) G R_ and |z| < 1 are now parameterised by the curves z + (t) 
and Z-(t) for t g] 0,t am ] where t am is given by setting s = am in (3.6); we write this 
out explicitly for future reference: 



tam 



am/2 (3.8) 



\ am/4 + ^1 + (am/4) 2 

These curves, which we denote by C + and C_ , lie in the lower- and upper half-planes 
of C, respectively. They have a common limit point at z+(0) = Z-(0) = 0. See Fig.l. 

According to the residue theorem, the integral (3.2) remains unchanged when the 
unit circle is shrunk to a closed curve C around C + U {0} UC_. In the limit this 
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reduces to an integral over C + U {0} U C_ itself, with a factor of 2 to take account 
of the fact that C goes along C + U C_ twice, with opposite orientations. (This is 
assuming that the argument in z [X ~ v ^ a is sufficiently large to avoid a divergence of 
the limit integral due to singularity at z = ; an explicit criterion for this will be 
given further below.) 5 Then, using obvious symmetries, the integral can be seen to 
be 2 x the real part of the integral over C_. (The integrals over C + and C_ are 
complex conjugate, so the imaginary parts cancel out as they should.) The square 



root y f(z) in the integral then reduces to ±iyX(t) with the explicit expression for 
X(t) determined below and the sign determined to be "— ". Thus (3.2) reduces to 



am) 2 + A 2 (x, y) = — 

na Jo 



dt 



dz- 



dt 



X(t) Refz-it^-M*- 1 ) (3.9) 



Recalling that the solution to f(z) = —X can be written as (3.5) with s = ^ (am) 2 + A 
and noting that the relation between s and t in (3.6) can be inverted to give s = 
2(1 - t) 2 1 y/t(2 - t) , we find 



4(1 - ty 

\l t(2 - 1) 



— (am)^ 



(3.10) 



The sign in ±iyX(t) can be determined by considering y f(z) ~ V z~ 2 for z near zero. 
Writing z = e + iS we have z~ 2 m (—5 2 — i2Se)/(5 4 + (2<5e) 2 ) ; hence, recalling that 
we have chosen the cut R_ to define the square root, z~ 2 —i/S. From this it is 
straightforward to see that the sign in ±i^JX(t) in the integral over C_ is "— ", and 
this is the origin of the minus sign in (3.9). Explicit expressions for the remaining 
ingredients in the integrand in (3.9) are readily found from (3.7): 



dz± 



dt 



l + t(2 -t) 
2-t \\ t(2-t) 



(3.11) 



(3.12) 



5 If the argument in z^ x ^ y ^ a is not sufficiently large, e.g. if x = y, then the curve around which 
the integral is performed cannot be completely shrunk to C + U {0} UC_ - a small detour around 
z = must be included. This case is more subtle, and we do not consider it here since it is not 
needed to derive the advertised locality result. 
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Note that the divergences ~ 1/ \fl for t — > in \J\(t) and ^ are compensated in 
(3.9) by powers of \ft in z_(£)l x-s/ l/° -1 provided \x — y\ > a, which we henceforth 
assume to be the case. (This is the criterion alluded to above.) We can now use (3.9) 
to draw conclusions about the exponential decay of yj (am) 2 + A 2 (x,y). Explicit 
evaluation of the integral in (3.9) will not be needed for this, so we do not attempt 
to perform it here. 

For fixed \x — y\ > and given t G [0, t am ] the integrand in (3.9) is dominated in 
the a — > limit by the exponential factor z-(t)^ x ~ y ^ a . From (3.11) we see that 
increases with t for t e [0, t am ] (recall t am < 1); therefore there can be no cancellation 
between the exponential factors for different t in the integral (3.9) and it follows that 



(am) 2 + A 2 (x,y) decays exponentially ~ z^(t am y x y ^ a for small a. From (3.7) we 
see that 



t am = 1 _ J am/2 + O(am) . (3.13) 



Consequently, using (3.11), the magnitude of the exponential decay of y (am) 2 + A 2 (x, y) 
for small a (i.e. am « 1) is found to be 

\z-(t am )\ lx - yl/a = (l-^/2 + 0(am)) M/2S 

iA/5 1 \ x ~y\/ 2 ^ 



(l - ym/2v^ + 0( am) 
%° (e-^ 2 f- yy2VE = e -i>/fi \*~v\ (3.14) 



Thus the localisation range for the exponential decay of y (am) 2 + A 2 (x, y) is seen 
to be 2y / 2a/m . 



We now supplement the preceding with a bound on |y (am) 2 + A 2 (x,y)\ which 
allows to check that the integral in (3.9) does not give rise to other factors which 
mask the exponential decay when \x — y\ is of the same order of magnitude as yja/m. 
From (3.10)-(3.11) we see that for t E [0,1] 



dz- 



dt 



< ^ , < , !*-(*) I < ^ (3-15) 
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and it follows from (3.9) that 

/ 4\/2 r ta ™ / r- \ 

J (am) 2 + A 2 (x,y) < — / dt(Vt) a 
v ira Jo v ' 



Sv 7 ? /. \\*-y\/^ , . 

1 Xitam) (3.16) 

TTV^am (F - J/| - aj 



For am « 1 the exponential factor here reduces as in (3.14) to give the same decay 
found earlier. The factor 1 / \Jt am has no effect since by (3.13) it is I. On the other 
hand, the factor 1/ (\x — y\ —a) blows up for \x — y\ ~ a; however it does not mask the 
exponential decay once \x — y\ > 2a (and enhances the locality when \x — y\ is large). 6 



When a is sufficiently small, the localisation range (~ ya/m ) of the exponential 
decay is much larger than a and therefore does not get masked by this factor. 

We now proceed to the case of arbitrary spacetime dimension d and consider 

JH 2 + EA^ (x,y) = -^— d [ d d p I (am)* + £ A u ( Pl/ y e'E.M^-v,)/. 
V v (2na) a J[-ir,ir]d Y - 

(3.17) 

Writing x = (xi, x) , p — (pi, p) and setting 

M(p) 



\ 



am) 2 + ^A^) 2 (3.18) 



^=2 



we have 



^(am) 2 + ]TA 2 (x,y) 

= - — 1 —— ( d d ~ l v e ip ' x/a r JM(p) 2 + A t ( Pl ) 2 e iMxi - yi)/a 

(27ra) rf - 1 J[-«,n]d-i V J-n2na y 1KF1J 

(3.19) 

The integral over pi here is the same as the previous d — 1 integral (3.2) except that m 
is replaced here by M(p). It can therefore be rewritten as (3.9) with this replacement. 
By our previous argument this integral decays exponentially ~ 2-(£M(p))' Xl-s/1 '^ a - The 
decay is slowest when £m( p ) is largest, i.e. when M(p) is smallest, and this happens 
when p = (0, . . . , 0) in which case M = M(0) = am. The same reasoning which led 
6 Recall that the derivation of (3.9), and hence also (3.16), assumes \x — y\ > a. 
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to (3.14) then implies that for am « 1 the operator kernel ^ {am) 2 + J2u { x iV) 
decays ~ e~2\/i^ l^ 1- ^ 1 ! along the — 1 axis. Obvious modifications in the preceding 
show that the same decay holds along any other coordinate axis. Thus we see that 
the localisation range is no smaller than 2v / 2y / a/m. It could however be larger along 
directions which are not parallel to a coordinate axis. To derive an upper bound on 
the localisation range we use bounds similar to those leading to (3.16) to get 7 

l/^f^M ^s^fc^W'"' 1 ''' (3 - 20) 

holding for each /i — 1, 2, . . . , d. It follows that 

L/w+Eai (»,»)f < n *-*f — nW"""* 

V v "i a ^Am (F/i-J^I -a) V > 
which in turn gives 

L/h'+eam^)! < ( tm f-*>™ 

V " a d ~H^U^ (n^d^-y^l -a) J 

(3.21) 

where — y\\\ = ^« k// ~ is the 'taxi-driver' norm. A calculation analogous to 
(3.14) gives 



(f am ) l|x " wlll/MB a 3° e "Wfi H*"^ . ( 3. 22) 
Since ||x — y\\ < ||x — y||i it follows from this and (3.21) that the localisation range 



is no bigger than 2d^2a/m, i.e. it lies between this value and the previously derived 
lower limit 2^2a/m. This completes the demonstration of exponential locality, with 
localisation range vanishing ~ yja/m , claimed at the beginning of this section. 

It is interesting to compare this result with the free field locality result derived in 
[9] for the operator 



(am) 2 + ^(V,)tV,. (3.23) 



7 The factor originates from the first integral in (3.19): ^aV" 1 n]*- 1 d d 1 P |e 4p x / a | 
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This operator was shown there to be exponentially local but with localisation range 
remaining finite in the a — > limit. The argument involved a mixture of analytic and 
numerical calculations 8 ; however, the result can be established by purely analytic 
means, using the techniques introduced in the preceding, as we now demonstrate. 
This will also show the origin of the difference in locality behaviour between our 
operator and this one. (Note J2 U Vj, = — V 2 ; we use the latter expression in the 
following.) 

In the d= 1 case, 

J (am) 2 - V 2 (x, y) = T dp J (am) 2 - V 2 (p) e tp{x - y)/a 

= t^— I — J (am) 2 -(z- z' 1 ) 2 z\ x ~ y \ /a (3.24) 

Setting g(z) = (am) 2 — (z — z~ r ) 2 = z~ 2 {(am) 2 z 2 — (z 2 — l) 2 ^) we proceed as before 
by determining the z's satisfying g(z) G R_ and \z\ < 1 : 

g(z) = -A A e R+ 
& -(z 2 -If + ((amf + \)z 2 = 

<£> (z 1 - 1 + ^j(am) 2 + A z) (z 2 - 1 - ^ (am) 2 + A 2) = (3.25) 

The solutions with \z\ < 1 are z± (s = ^(am) 2 + A ) where 

z+(s) = -s/2 + y/l + (s/2) 2 = -z-(s) s E [0,oo[ (3.26) 

Note that ^±(s) — >■ for s — > 00. Hence the solutions form curves C + and C_ inside 

the unit circle in C, parameterised, respectively, by z + (s) and Z-(s) , s e [am, 00 [. 

The curves in this case are simply intervals on the real axis: C + =] 0, z + (am)] and 

C_ = [-z + (am),0[ (see Fig. 2). The circle around which the integration in (3.24) is 

carried out can now be shrunk to a closed curve around C + U {0} U C_ , leading in the 

limit to an integral over these curves. By arguments similar to those in the previous 
8 Specifically, the locality behaviour of the continuum version of this operator was analytically 
determined and numerical calculations were then performed to check that the lattice operator kernel 
reduced to the continuum expression in the a — > limit - see Part 3 of Appendix B in [9]. 
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Figure 2: The "exclusion curves" C± for the operator (3.24). Compare with Fig. 1 
for the operator (3.2). 

case one then finds that yj (am) 2 — V 2 (x,y) decays exponentially ~ \z + (am) \\ x ~y\/ a . 
From (3.26) we see that for am « 1 the magnitude of the decay factor becomes 

\z + (am)\ lx - yl/a = ((l-am/2 + 0(H 2 )) 1 /j M 

~° e"¥M. (3.27) 

Thus the localisation range in this case, 2/m, is independent of a and remains finite in 
the a — + limit. The general dimension d case can now be dealt with by an argument 



analogous to our earlier one; this leads to the result that y (am) 2 — V 2 (x, y) is 
exponentially local with lower- and upper bounds on the localisation range being 
2/m and 2d/m, respectively, showing that the range is also finite ~ m _1 in the a — > 
limit in the general dimension d case. Thus we reproduce the general finding of [9] 
for the locality behaviour of this free field operator. 9 

The origin of the different locality behaviour of our operator and the free field 

operator considered in [9] is now apparent: The exponential decay in d— 1 dimension, 

which, as we have seen, is the same as the decay along a coordinate axis in general 
9 In the expression for the free field operator kernel G(x,y) in Eqn.(3.11) of [9] the integration 
range for dp^ (after a change of variables p^ — > p^/a) is [ — ir/2, n/2]. But since — V^(p M ) = 
sin 2 (p^) this gives precisely 1/2 of what the integration over [— 7r,7r] would give. Hence G(x,y) = 



■^yj(am) 2 — V 2 (x, y) so the locality result derived above applies. 
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d dimensions, is given in both cases by \z max \^ x ~ y ^ a , where z max is the point on the 
"exclusion curves" in Fig.l (our case) or Fig.2 (the case of Ref.[9]) which is closest to 
the unit circle. In our case, 

\z ma x\ ~ (t am ) 1/2 w (l- y/am/2 + O (am)) ^ , (3.28) 

resulting in decay ~ e'^V^ \ x ~ y \ , whereas in the case of Ref.[9], 

\zmax\ « 1 - am/ 2 + 0{{amf) , (3.29) 

resulting in decay ~ \ x ~ y \. The essential difference is that the leading a-dependent 
term inside the bracket in (3.28) is ~ ^Jam whereas the corresponding leading term 
in (3.29) is ~ am. This leads to the localisation range being ~ ^ a/m in the former 
case and ~ m (independent of a) in the latter. The former tends to zero for a — > 
while the latter stays constant in this limit. 

We remark that the technique of writing the kernel 0(x,y) of a free field lattice 
operator as an integral around a closed curve in the complex plane, and then attempt- 
ing to shrink the curve as a way of deriving locality properties, was used previously 
in a different context in Ref. [24] . 

4 Discussion 

Working in the flavour (taste) representation for staggered fermions, we have shown 
that in the free field case there is a simple and natural Wilson-type lattice Dirac oper- 
ator D on the blocked lattice, given by (2.10), which satisfies detD = {detDgtagg^.^ 1 ^ 
and is exponentially local with localisation range vanishing ~ ^ a/m for a — > 0. The 
techniques developed to derive the free field locality result are of a generally appli- 
cable nature, and we also used them to give a new, purely analytic derivation of the 
negative locality result in [9]. They can also be used to study free field locality prop- 
erties of other lattice operators of current interest; in particular, the overlap Dirac 
operator [25] , which is treated in a forthcoming paper [26] . 
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Our free field operator can be gauged, i.e. coupled to the link variables of the 
original lattice, in a variety of ways. The simplest way is to define link variables 
on the blocked lattice in terms of the link variables on the original lattice by 

V^2x) = U^2x)U^2x + afi) (4.1) 

(/t=unit vector in the positive /z-direction); then the difference operators , A u on 
the blocked lattice in (2.10) can be coupled to in the usual way - this specifies the 
"minimal gauging" of our D. However, the resulting operator cannot be expected to 
satisfy the SDR in the interacting case. Our argument in §2 does not carry over to 
this case; it is specific to the free field case. 

Regarding the possibility of gauging our operator such that the SDR does continue 
to hold in the interacting case, we note the following. The taste-mixing part of the 
staggered Dirac operator leaves unbroken a U(l) subgroup of the continuum U(4) axial 
flavour symmetry. 10 In the taste representation, this symmetry can be expressed in 
the free field case as 

{ l5 ®T 5 ,D stt } = (m = 0) (4.2) 

with notations as in (2.2). This chiral symmetry protects staggered fermions against 
additive mass renormalisation [27, 28, 29]. The sea quark effective action in dynamical 
staggered fermion simulations is log det{D st ) 1 ^ = ^Tr logD st , the same as for a usual 
staggered fermion modulo an overall factor 1/4. Thus the protection against additive 
mass renormalisation is also present in this case. This situation would be difficult, 
if not impossible, to reconcile with the existence of a single-flavour D satisfying the 
SDR unless the D is also protected against additive mass renormalisation: if the 
bare mass is small then for staggered fermions the physical mass will also be small, 
whereas for a lattice fermion described by a D which does not have a chiral symmetry 
10 We are assuming that the mass term of the staggered fermion is of the standard form m(l ® 1). 
(If the mass matrix is not proportional to the identity operator then the interpretation of the U(l) 
symmetry is different.) Note that the U(l) symmetry is not the diagonal U(l) subgroup in U(4). The 
latter, associated with the axial anomaly, is explicitly broken by the taste-mixing part of D sta ggered- 
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the physical mass will be large due to the additive mass renormalisation induced by 
radiative corrections. Since our D is of Wilson-Dirac form, any gauged version of 
it will be afflicted with additive mass renormalisation, unless there is a very special 
choice of gauging which endows this D with a new, hitherto undiscovered type of 
chiral symmetry. The latter seems very unlikely though, so most probably a gauged 
version of our D satisfying the SDR simply does not exist. 

The symmetry corresponding to (4.2) for (D stt )^D stt in the free field case is 

[ 7 5 <g) T 5 , (D stt yD stt \ = (m = 0) . (4.3) 

But, as noted in §2, (D stt yD stt ~ 1 <g> 1 in the free field case, so (4.3) is trivially 
satisfied. This explains why it was possible to find a single-flavour D without chiral 
symmetry but nevertheless satisfying the SDR in the free field case. In the interacting 
case (D stt y D stt is no longer ~ 1 <g> 1 and the gauged version of the symmetry (4.3) is 
a nontrivial property. 

Although a gauged version of our D satisfying the SDR is unlikely to exist, the 
free field locality result for it is still relevant as a general indication of the possibility 
of having local single-flavour Dirac operators satisfying the SDR, and as a first step 
toward constructing more sophisticated operators which have a chance to be local 
(also for m = 0) and satisfy the SDR in the full interacting case. At present it is the 
only analytic positive locality result derived for a solution of the SDR in any gauge 
background (the background in our case being the trivial one). 

The above discussion indicates that, for a viable D to satisfy the SDR in the full 
interacting case, it should have an exact chiral-type symmetry (presumably corre- 
sponding in some way to the aforementioned chiral symmetry of staggered fermions). 
The only such symmetry currently known for single-flavour lattice Dirac operators 
is the lattice-deformed chiral symmetry [30] possessed by operators satisfying the 
Ginsparg- Wilson (GW) relation [31, 32, 33] and its generalisations [34]. This sug- 
gests to look for free field solutions to the SDR which also satisfy the GW relation, in 
the hope that among these there may be a D which can be gauged such that the SDR 
continues to hold in the interacting case. In fact this was already investigated by 
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numerical means on finite lattices in Ref.[13]. The numerical results there appear to 
be encouraging. However, the problem can also be addressed analytically and when 
this is done difficulties become apparent. Setting the lattice spacing of the blocked 
lattice to unity for convenience, the GW relation in its broad form is 

75D + £75 = 2D lb RD (4.4) 

where R is an arbitrary local scalar hermitian operator. As noted at the end of §2, a 
sufficient condition for the free field D to satisfy the SDR is 

D^D = A + m 2 (4.5) 

In the free field case, assuming 75-hermiticity = 75-D75 and exploiting the fact 
that solutions to (4.4) are of the form D = (2R)~ 1 (1 + 7 5 e) where e 2 = 1 and e' = e , 
the most general solution to (4.4) is seen to be of the form 



D. 



GW 



(4.6) 

m v v -v 2 + w 2 J 



where and W are sums of scalar operators multiplied by an even number of 7- 
matrices. Straightforward algebra now shows that requiring this operator to satisfy 
(4.5) fixes the W such that 



D G w = 7^^ 



(A + m 2 )(l- ±(2i?) 2 (A + m 2 )) 



+ R(A + m 2 ) (4.7) 



\R\\ ~ -V 2 

This is the general solution to (4.4)-(4.5) in the free field case. The numerical solu- 
tions investigated in Ref. [13] are particular cases of this operator (or more precisely, 
approximations to it on finite lattices) with the V M 's being scalar operators. The 
main interest here is in the case m — 0, since it is in the chiral limit that D should 
have the chiral symmetry implied by the GW relation (4.4). 

In the simplest case where R = 1/2 and = , i.e. the usual symmetrised 
difference operator, D GW reduces in the m = case to 



1 - \ E (^) + |A (4.8) 
20 



While this operator correctly reproduces the continuum free field Dirac operator in the 
a — > limit, it is most unlikely to have acceptable locality behaviour. The presence 
of the (— V 2 ) -1 inside the square root in the general operator (4.7) makes it difficult 
to envisage that there exist V M 's and a local R for which this operator has acceptable 
locality behaviour either, in spite of the numerical indications from Ref. [13]. Thus 
it would seem that the condition (4.5), which is sufficient, but not necessary, for the 
SDR to be satisfied, is actually too restrictive to lead to a local operator D satisfying 
both the SDR and GW relation. 

The preceding considerations indicate that finding a viable exact solution to the 
SDR in the interacting case is a difficult problem. However, to resolve the fourth 
root issue is is not actually necessary to have an exact solution; it suffices to find a 
viable lattice Dirac operator which satisfies the SDR approximately in the sense that 
the effective action difference 

d(U) = \ log detD st - log detD (4.9) 

is effectively just a lattice Yang-Mills action for the gauge field. In this case, represent- 
ing the quark determinant by det(D st ) 1 ^ is physically equivalent to representing it by 
detD together with a renormalisation of the bare coupling constant (i.e. a shift in (5). 
In other words, d(U) has the same effect as the fermion determinant for dynamical 
heavy quarks [21]. In connection with this it is useful to note that the perturbative 
expansion of a general single-flavour lattice fermion determinant has the form [35] 

oo 

log detD = (—^z log(am) 2 + c D ) S Y m{A) + ^2(I n (A;m) + v n (A; am)) 

n=2 

(4.10) 

where Sym(A) is the continuum Yang-Mills action, I n (A; m) is a non-local continuum 
functional of order n in A, and the v n (A; am)'s (also nonlocal and of order n in A) 
are terms which vanish for am — > 0. The dependence on the specific choice of lattice 
Dirac operator D enters only through the numerical coefficient cp and the functions 
v n (A;am). (A gauge field-independent term which diverges for am — > has been 
ignored in (4.10).) In fact the perturbative expansion of \log detD st has the same 

21 



form (4.10) as a single-flavour Dirac operator [35]. Letting c st denote the coefficient c D 
in this case, it follows that the perturbative expansion of the effective action difference 
(4.9) has the form 

oo 

d(A) = (c st - c d )Sym(A) + w n (A; am) (4.11) 

n=2 

with each w n (A; am) vanishing for am — > 0. Thus it would seem that, in the pertur- 
bative setting at least, for am « 1 the effective action difference is indeed just a 
Yang-Mills action for the gauge field, for any sensible choice of single-flavour D. The 
situation is not this simple though - although they vanish for a — > 0, the functions 
v n in (4.10) and w n in (4.11) still do affect the quantum continuum limit. At the 
perturbative level this is manifested in that Feynman diagrams with vertices from 
these terms can be non- vanishing; in fact divergent. To see this, recall that the 
terms in the perturbative expansion of log detD are given by l-(fermion)-loop gluonic 
n-point functions. The internal propagators and vertices in these receive radiative 
corrections. In particular, unless D is protected by a chiral symmetry, the radiative 
corrections to the internal propagators give rise to a large additive mass renormalisa- 
tion. Since there is no corresponding effect from D st to cancel this, it will manifest 
itself in the divergence (for am — > 0) of various Feynman diagrams involving vertices 
from the "irrelevant" terms w n (A;am) in the effective action difference (4.11). Thus 
the importance of D having a chiral symmetry becomes clear in this context as well. 

In perturbative (lattice) QCD the renormalisations of interaction vertices are inde- 
pendent of the renormalisation of the fermion propagator and bare mass. Thus (4.11) 
and the observations above would suggest that, for am « 1, the effective action dif- 
ference d(A) is indeed essentially a Yang-Mills action when the single-flavour lattice 
Dirac operator D has an exact chiral symmetry; e.g., when D is the overlap Dirac 
operator. Support for this hypothesis comes from a numerical study carried out in 
two dimensions in Ref. [17]. In two spacetime dimensions the perturbative expansion 
of log detD is completely universal, modulo terms which vanish for am — > [35] (this 
is a reflection of the fact that QCD in two dimensions is super-renormalisable). Thus 
the first term in the right-hand side of the effective action difference (4.11) is absent 
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in this case, and the hypothesis then states that representing the quark determinant 
by det(D st ) l l 2 is equivalent to using detD ov (with D ov being the overlap Dirac opera- 
tor) without any renormalisation of the bare coupling. If this hypothesis holds, then 
det(Dst) 1 / 2 should coincide with detD ov for equilibrium gauge configurations of an en- 
semble generated by taking the probability weight to be e~^ SyM ^ detest) 1 / 2 . And 
this is precisely what was found to good accuracy in a numerical study in Ref. [17]. 

It must be remembered though that the perturbative picture is not the full picture. 
Low-lying eigenvalues of the Dirac operator are associated with long-range, low energy 
dynamics in QCD which is not captured by the perturbative framework. Indeed, 
numerical studies of the Wilson fermion determinant in Ref. [36] show that the log of 
the determinant cannot be modelled by a linear combination of local loop functionals 
(i.e. functional of the form Tr U(a) where U (a) is the product of the link variables 
around a closed lattice path a) ; in particular it cannot be modelled by a local lattice 
YM action. However, the product of the Dirac eigenvalues of magnitude > Aqcd does 
admit such a description [36] . Thus, the aforementioned hypothesis, coming from the 
perturbative considerations above, should be regarded as applying to the truncations 
of det{D st ) l / A , detD, and the effective action difference d(U), given by excluding the 
eigenvalues of magnitude < Aqcd- (A way to implement and study this truncation 
in the perturbative setting is mentioned in [37].) 

Specifically, defining detD high and detD[ ow to be the products of the eigenvalues 
of D of magnitudes > Aqcd, and < Aqcd, respectively, and splitting up the effec- 
tive action difference into d = dhigh + di ow in the obvious way, the hypothesis can be 
stated as follows: "When D is the overlap Dirac operator and am« 1 then d high (U) 
is essentially a local lattice YM action in 4 dimensions, and essentially vanishing in 
2 dimensions." The question of whether d(U) itself is effectively a local YM action 
(which would give a positive resolution of the fourth root issue if the answer is affir- 
mative) is then reduced to a question of whether or not di ow {U) is effectively zero. 
Thus it would be highly desirable to numerically study det{D st ^ ow ) 1 ^ , detD ov j ow , and 
thereby di ow (U), in equilibrium gauge backgrounds in 4 dimensions. We remark that, 
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in 2 dimensions, combining the hypothesis that d high vanishes (in 2 dim.) with the 
numerical agreement [17] between the full rooted staggered and overlap determinants 
implies that di ow does indeed vanish in this case. 

In 4 dimensions numerical studies have found that, after applying a UV- filtering 
procedure, there is good agreement between the low-lying eigenvalues of D ov and D st 
(modulo a four-fold degeneracy in the latter) [17]. 11 Comparisons of the spectrum of 
D st with predictions of Random Matrix Theory also back up this picture [18, 19]. 
While this does not by itself prove that di ow vanishes, it is certainly compatible and 
suggestive of it. 

While the numerical work in this direction may lead to a resolution of the fourth 
root issue at a practical level, one should ask whether it is possible to also get a 
resolution at the theoretical level in this approach. If it is possible it will probably 
happen as follows: (i) Use renormalisation group arguments to justify a perturbative 
treatment of dhigh and verify the hypothesis that it is effectively a local lattice YM 
action, (ii) By applying Random Matrix Theory and theoretical implications of UV- 
filtering to the low-lying spectra of D ov and D st show that d\ ow is effectively zero when 
taking the quantum continuum limit. 

Another interesting and promising approach to the fourth root issue has been 
given recently by Shamir [14]. A renormalisation group argument is used to express 
the free field staggered fermion action in the flavour (taste) representation on a lattice 
spacing a as an action on a coarse lattice of spacing a = 2 n ao . This results in a 
decomposition of the staggered fermion determinant in the form 

det(D st ) = det(D n ) de^G' 1 ) (4.12) 
11 In fact it is only when this UV-filtering is applied that the aforementioned agreement between 
the rooted staggered and overlap determinants in 2 dim. holds [17]. Without the filtering the 
agreement breaks down, just as the agreement between the low-lying eigenvalues does. So it is 
tempting to ascribe the breakdown in the agreement between the determinants to the breakdown in 
the agreement between the low-lying eigenvalues. This gives a further hint that vanishing of d low is 
intimately connected with having agreement between the low- lying eigenvalues of D ov and D st . 
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The operator D n encodes the low energy/long range dynamics of staggered fermions; 
it decays exponentially with localisation range ~ a , satisfies a GW relation when 
m — 0, and becomes proportional to the identity matrix in flavour space in the large 
n limit: lim^oo D n = D rg <g> 1 flavour ■ Hence, in this limit 

detiDst) 1 ^ = det(D rg ) det{G^ 1/4 ) . (4.13) 

The magnitude of the spectrum of each G^ 1 has a lower bound ~ a; in units of the 
fine lattice spacing a this blows up for large n, so the expectation is that in a gauged 
version of this setting the effect of det{G^/ A ) in (4.13) is exactly the same as that 
of the determinant of a heavy dynamical fermion: to simply renormalise the bare 
coupling parameter. 

It should be pointed out though that the n — > oo limit leading to (4.13) cannot 
actually be taken in practise - it corresponds to a — > 0, but a must remain non-zero 
since it is the spacing of the lattice on which the staggered fermion lives and the 
Lattice QCD simulations are performed. Therefore, in this approach one needs to 
remain at finite n, i.e. the setting of (4.12). For large finite n the operator D n is close 
to being diagonal in flavour space, but is not exactly diagonal. This is different from 
the situation in the present paper where we obtain a single-flavor candidate Dirac 
operator already at non-zero lattice spacing. To fully resolve the fourth root issue 
in Shamir's approach it is necessary to find an single-flavour lattice Dirac operator 
D' such that adding log detD' — \log detD n to the lattice gauge field action does not 
affect the quantum continuum limit. Shamir has a proposal for this operator D' [38]. 
Moreover, his approach has a definite possibility of being extended to the interacting 
case, although this remains a difficult challenge for future work. 

An appealing feature of Shamir's approach is that the GW chiral symmetries of D n 
and D Tg (at m — 0) originate in a clear and direct way from the chiral symmetry (4.2) 
of the staggered Dirac operator. This also raises intriguing questions. The chiral 
symmetry of a GW Dirac operator is generally anomalous - it gets broken by the 
fermion integration measure [30] . On the other hand, the chiral symmetry of D st gets 
broken spontaneously in the m — > limit (at least at strong coupling), and there is an 
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associated Goldstone meson [39, 40]. In connection with this we mention a potentially 
troubling aspect of the fourth root prescription which has been pointed out already by 
Creutz [41]: When the determinant for a single quark is represented by det(D st ) 1 ^ 4: , 
what becomes of the Goldstone meson associated with the spontaneous breaking of 
the chiral symmetry of D st l Single-flavour Dirac operators are not supposed to have 
spontaneously broken chiral symmetries. This and other intriguing issues for the 
fourth root prescription remain as an urgent topic for future work. 

Finally, we mention that a completely different approach to this issue, involving 
relating the fourth root prescription to local theories via a parameter deformation in a 
family of lattice theories in 6 spacetime dimensions, has been described by Neuberger 
in Ref.[12]. 
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